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Computation  of  Forces  on  a  Rotating  WindSat  Reflector 


INTRODUCTION 

A  satellite  based  wind  speed  and  direction  detecting  system  called  the  WindSat  has  been 
developed  by  NRL,  The  main  reflector  of  the  WindSat  system  rotates  and  is  connected  to  the 
Coriolis  spacecraft  via  an  interface  and  stationary  platform.  Any  unbalanced  force  produced  by 
the  rotating  antenna  would  result  in  a  wobble  in  the  entire  spacecraft.  It  is  therefore  critical  to 
assess  the  forces  produced  by  the  rotating  antenna  and  counterbalance  appropriately.  In  order  to 
measure  imbalance  of  forces,  the  reflector  assembly  is  tested  with  the  test  bed  spinning  in  the 
clockwise  direction.  The  direction  of  rotation  is  then  reversed  and  the  forces  are  measured. 
From  these  two  tests,  the  imbalance  of  forces  in  the  circumferential  direction  and  its  moments 
about  the  spin  axis  can  be  obtained  and  corrected  by  adding  counterbalancing  weights.  The 
force  normal  to  the  reflector,  a  lift  force,  and  its  moments  about  the  axis  of  rotation  can  not  be 
determined  from  these  two  tests.  The  aim  of  this  effort  is  to  computationally  obtain  the  lift  force 
on  the  rotating  reflector. 


INCOMPRESSIBLE  FLOW  SOLVER 

FEFLO  is  a  finite  element-based  incompressible  flow  solver  based  on  simple,  low-order 
elements.  The  simple  elements  enable  the  flow  solver  to  be  as  fast  as  possible  reducing  the 
overhead  in  building  element  matrices,  residual  vectors  etc.  The  governing  equations  are  written 
in  Arbitrary  Lagrangian  Eulerian  form  which  enables  simulation  of  flow  with  moving  bodies. 
For  high  Reynolds  number  flow  cases,  the  mesh  requirement  is  met  by  employing  arbitrary  semi- 
structured  grids  close  to  wetted  surfaces  and  wakes.  The  full  details  of  the  flow  solver,  the  rigid 
body  motion  and  adaptive  remeshing  are  given  by  Ramamurti  et  al.1  and  are  summarized  next. 

The  governing  equations  employed  are  the  incompressible  Navi er- Stokes  equations  in 
Arbitrary  Lagrangian  Eulerian  (ALE)  formulation.  They  are  written  as 


Dv 

—  +  va  *Vv  +  Vp  =  V*tr  (la) 

Dv  dv 

—  -  —  +w-Vv  (lb) 

Dt  dt  ' 

V  ■  v  =  0  (2) 

Here  p  denotes  the  pressure,  va  =  v-  w,  the  advective  velocity  vector  (flow  velocity  v  minus 
mesh  velocity  w),  and  the  material  derivative  is  with  respect  to  the  mesh  velocity  w.  Both  the 
pressure  p  and  the  stress  tensor  a  have  been  normalized  by  the  (constant)  density  p,  and  are 
discretized  in  time  using  an  implicit  time  stepping  procedure.  It  is  important  for  the  flow  solver 
to  be  able  to  capture  the  unsteadiness  of  a  flow  field.  The  present  flow  solver  is  built  as  time- 
accurate  from  the  onset,  allowing  local  timestepping  as  an  option.  The  resulting  expressions  are 
subsequently  discretized  in  space  using  a  Galerkin  procedure  with  linear  tetrahedral  elements.  In 
order  to  be  as  fast  as  possible,  the  overhead  in  building  element  matrices,  residual  vectors,  etc. 
should  be  kept  to  a  minimum.  This  requirement  is  met  by  employing  simple,  low-order  elements 
that  have  all  the  variables  (u,v,w,p)  at  the  same  location.  The  resulting  matrix  systems  are  solved 
iteratively  using  a  preconditioned  gradient  algorithm  (PCG).  The  flow  solver  has  been 
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successfully  evaluated  for  both  2-D  and  3-D,  laminar  and  turbulent  flow  problems  by  Ramamurti 
etalP 

In  a  frame  of  reference  rotating  with  an  angular  velocity  Q,  Eq.  (la)  can  be  written  as 

Oy 

—  +  vfl  *Vy  +  Vp  -  V-  a-2Qx  v-Qx(Qx  v)  (3) 

Compared  to  the  stationary  frame  of  reference,  we  have  additional  source  terms  from  the 
Coriolis  and  centrifugal  forces. 

COMPUTATIONAL  RESULTS 

The  FEFLO  flow  solver  in  the  rotating  frame  of  reference  was  first  validated  using  two  test 
problems.  The  first  test  problem  is  the  flow  near  a  rotating  disk.  The  flow  in  the  vicinity  of  the 
disk  is  validated  with  exact  solution.  This  is  essential  for  computing  forces  and  moments  on  a 
reflector-like  geometry  where  the  pressure  on  the  surface  and  the  velocity  in  the  vicinity  of  the 
disk  are  the  most  important  parameters.  The  second  test  case  is  that  of  flow  between  two  rotating 
cylinders.  The  purpose  of  this  test  is  to  validate  the  computations  in  the  rotating  frame  of 
reference. 


Flow  near  a  rotating  disk 

The  first  test  problem  is  the  flow  near  a  rotating  disk  shown  in  Fig.  1.  A  flat  plate  rotates 
about  an  axis  perpendicular  to  its  plane  with  a  uniform  angular  velocity  w,  in  a  fluid  otherwise  at 
rest.  An  exact  solution  of  the  Navier-Stokes  equations  is  possible  for  this  case  and  is  given  by 


Schlichting4.  Introducing  a  dimensionless  distance  from  the  wall  £  =  and  using  the 


following  assumptions  for  the  radial,  circumferential  and  axial  velocity  components  and 
pressure, 

u  =  ra)  F(Q,  v  =>  roi  G(£),  w  =  -y[va)H(£)  and 
p  =  p(z)  =  pvw  P(Q  , 

Eq.  (1)  can  be  written  as  a  system  of  four  simultaneous  ordinary  differential  equations  for  the 
functions  F,  G,  H  and  P: 

2F  +  /T  =  0 


F2  +  F'H-G2-F*  =  0 
2FG  +  HG'  -  G*  =  0 


(5a) 


P'  +  HH'-H"  =  0  . 

The  boundary  conditions  for  this  system  can  be  obtained  from  the  no  slip  conditions  at  the  wall 
and  no  flow  far  away  from  the  plate  and  are: 

at  £  =  0 :  F=0,  G=l,  H=0,P=0  and 

(5b) 

at  %  =°o;  F=0,  G=0. 

Equation  5  was  solved  numerically  and  the  velocity  components  were  then  computed  using 
Eq.(4). 


Computations  of  Forces  on  WindSai  Reflector 


3 


The  FEFLO  flow  solver  was  then  used  to  solve  the  flow  past  this  configuration  as  shown  in 
Fig.  2.  Here  the  flat  plate  was  of  a  finite  diameter  of  1  and  the  far-field  boundaries  were  set  to  be 
at  5.  The  computed  velocities  along  a  line  at  z-  0.4,  shown  by  the  dotted  line  in  Fig  2,  were 
compared  to  the  exact  solution  of  Eq.  (5),  and  are  shown  in  Fig.  3.  There  is  excellent  agreement 
with  the  theory  close  to  the  wall.  The  radial  and  axial  components,  u  and  w  respectively,  are  in 
exact  agreement.  There  is  a  small  discrepancy  in  the  tangential  component  (v)  away  from  the 
flat  plate.  This  is  due  to  the  nearness  of  the  outer  boundaries  and  also  due  to  the  fact  that  the 
spinning  disk  is  of  finite  radius  and  the  edge  effect  cannot  be  neglected.  Figure  4a-d  show  the 
distribution  of  the  velocity  components  and  the  pressure  along  xz  plane  passing  through  the 
center  of  the  disk.  The  radial  («)  and  circumferential  (v)  velocity  components  distributions  are 
anti-symmetric  about  the  y  axis,  Fig.  4a  and  b.  The  effect  of  the  edge  and  the  outer  boundary  can 
be  seen  in  these  contours  near  the  edge  of  the  disk  for  small  values  of  z  and  the  effect  spreads 
away  from  the  wall.  The  effect  of  the  outer  boundary  is  most  pronounced  in  the  axial  (w) 
component  of  the  velocity.  Fig  4c.  The  pressure  distribution  shown  in  Fig.  4d  depends  on  the 
distance  from  the  wall  for  very  small  values  of  z,  consistent  with  the  boundary  layer  theory. 
However,  away  from  the  wall  a  radial  pressure  gradient  is  present  which  is  produced  to  the  effect 
of  the  edge. 


Flow  between  rotating  cylinders 

In  order  to  validate  the  FEFLO  flow  solver  in  rotating  frame  of  reference,  the  flow  between 
two  concentric  rotating  cylinders  was  chosen.  An  exact  solution  of  the  Navier-Stokes  equation  is 
given  by  Schlichting4.  In  this  example,  the  two  concentric  cylinders  move  at  different  but  steady 
rotational  speeds.  The  inner  and  outer  radii  are  r,  and  r2,  and  move  with  angular  velocities  to, 
and  coj,  respectively.  The  circumferential  velocity  v  is  given  by 


v(r) 


r2  -n  l 


r(w2r/  -  corf)  -  -  o>,) 


(6a) 


For  the  special  case  when  the  inner  cylinder  is  rotating  and  the  outer  one  is  at  rest,  the  velocity  is 


v(r)  =  (o. 


(6b) 


The  configuration  for  this  computation  is  shown  in  Fig.  5a  with  the  inner  radius  of  1  and  outer 
radius  of  3.  Computations  were  performed  in  both  stationary  and  rotating  frames  of  reference 
with  the  rotational  velocity  of  the  inner  cylinder  w  -  60  rpm.  Figure  6  shows  the  comparison  of 
the  circumferential  velocity  between  the  cylinders  for  both  of  these  computations  with  the  exact 
solution  (Eq.  6b),  and  the  agreement  is  excellent.  It  was  found  that  the  solution  in  the  rotating 
frame  of  reference  converged  very  slowly  mainly  due  to  the  explicit  treatment  of  Coriolis  force 
terms.  Figures  7  a  show  the  circumferential  velocity  distribution  on  an  xz  plane  and  is  anti¬ 
symmetric  about  the  origin.  The  pressure  and  the  magnitude  of  velocity  on  a  z  =  constant  (=0.5) 
plane  is  shown  in  Figs.  7b  and  7c,  respectively  and  are  ax i symmetric. 


WindSat  Reflector 

Having  successfully  validated  the  flow  solver  for  the  two  cases,  the  FEFLO  was  applied  to 
compute  the  flow  past  the  WindSat  reflector.  First,  the  reflector  was  considered  without  the 
supporting  trusses.  The  reason  being,  that  the  trusses  would  decelerate  the  flow  and  hence,  the 
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case  without  the  trusses  would  yield  a  conservative  estimate  of  the  lift  and  drag  forces  generated 
on  the  reflector.  The  geometry  of  the  WindSat  reflector  was  obtained  from  Code  8230  and  is 
shown  in  Fig.  8.  The  reflector  is  obtained  from  the  intersection  of  a  parabolic  surface,  x3  =  4  Cz, 
with  a  cylinder  of  71 .6  in  diameter.  The  axis  of  the  cylinder  is  offset  from  the  axis  of  the  parent 
parabola  by  50  in.  The  center  of  the  chord  in  the  plane  of  the  resultant  parabolic  section  is  then 
allowed  to  rotate  about  an  axis  offset  by  20in.  For  computations,  this  geometry  is  set  in  a 
rotating  frame  of  reference  with  the  outer  boundaries  at  15  ft  from  the  axis  of  rotation.  The 
boundary  conditions  on  this  far  field  are  prescribed  velocities  rotating  in  the  direction  opposite  to 
the  rotation  of  the  reflector. 

The  incompressible  flow  solver  described  above  was  employed  to  compute  the  steady 
inviscid  flow  about  this  configuration.  From  our  past  experience,  an  inviscid  flow  model  is 
sufficient  to  obtain  the  correct  lift  forces,  if  the  flow  over  this  configuration  is  not  separated.  The 
effect  of  viscosity  will  increase  the  drag  forces,  in  the  circumferential  direction.  The  effect  of  this 
drag  force  on  the  unbalanced  moments  can  be  obtained  by  reversing  the  rotational  direction.  The 
rotational  speed,  cu,  of  the  reflector  was  first  set  to  60  rpm.  The  steady  state  flow  over  this 
configuration  was  obtained  in  2000  time  steps.  The  force  on  the  reflector  surface  and  the 
moment  about  the  origin  were  obtained  by  integrating  the  pressure  distribution  on  its  surface. 
The  computed  force  components  for  this  case  are/*  =  -0.99  lbs,/,  =  3.398  lbs  and/  =  0.58  lbs. 
The  moments  about  the  local  origin  shown  in  Fig.  8  are  Mx  -  -6.01  ft-Ibs,  My  ~  -3.185  ft-lbs  and 
8. 112  ft-lbs. 

Figures  9a  and  b  show  the  pressure  distribution  on  the  top  and  bottom  surfaces  of  the 
reflector.  It  can  be  seen  that  regions  of  high  pressure  exist  on  the  edge  of  the  reflector  in  the  ~z 
direction  and  almost  half  of  the  bottom  surface  in  the  +z  direction.  Figures  10a  and  b  show  the 
distribution  of  the  magnitude  of  velocity  on  the  two  surfaces.  On  the  top  surface,  the  velocity 
distribution  is  almost  uniform,  but  on  the  bottom  surface  we  see  a  high  velocity  region  near  +x 
axis  or  the  leading  edge  and  a  small  recirculation  region  in  the  -x  direction  or  the  trailing  edge  of 
the  reflector.  Figure  1 1  shows  the  velocity  vectors  on  the  bottom  surface.  Again,  the  small 
recirculation  region  is  clearly  seen  in  Fig.  1  lb. 

Next,  the  flow  past  this  configuration  was  obtained  with  the  reflector  rotating  in  the  counter 
clockwise  direction  at  60  rpm.  The  force  and  moment  components  obtained  for  this  case  are 
/  =  -0.99  lbs,  /  =  3.382  lbs  and /  =  -0.58  lbs,  and  Mx  =  5.97  ft-lbs,  My  =  3.172  ft-lbs  and  Mz  = 
8.024  ft-lbs,  respectively.  This  shows  that  the  force  in  the  y  direction  remains  almost  the  same 
with  the  force  in  the  circumferential  direction  opposite  to  the  previous  case,  as  expected. 
Another  computation  was  performed  to  assess  the  forces  at  a  reduced  rotational  speed  of  «>  =  30 
rpm.  For  this  case,  the  computed  force  and  moment  components  are  /  =  -0.248  lbs,  /  —  0.849 
lbs  and  /  =  0.1637  lbs,  Mx  -  -1.51  ft-lbs,  My  =  0.8  ft-lbs  and  Mz  =  2.032  ft-lbs,  respectively. 
Comparing  the  magnitude  of  the  forces  to  the  case  of  co  =  60  rpm,  it  is  clear  that  the  force 
components  are  reduced  by  a  factor  of  4,  This  is  to  be  expected  as  the  forces  are  proportional  to 
the  square  of  the  velocity,  and  the  rotational  speed  was  decreased  by  a  factor  of  2. 

Effect  Of  Viscosity 

In  order  to  study  the  effect  of  viscosity  on  the  forces  developed  on  the  reflector,  the  viscous 
flow  over  the  reflector  was  computed.  The  Reynolds  number  (Re)  for  this  flow  based  on  the 
maximum  rotational  tip  velocity  and  the  diameter  of  the  reflector  is  1,2  x  IQ6.  A  turbulent  flow 
past  the  reflector  was  computed  using  a  Smagorinsky5  turbulence  model.  Figure  12  shows  the 
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convergence  history  for  the  forces  for  both  the  inviscid  and  viscous  turbulent  flow  simulations. 
The  converged  forces  and  moments  are  fx  =  -1,25  lbs ,  /  =  3.96  lbs  and  /  =  -0.664  lbs,  and 
Mx  =  -6.199  ft-lbs,  My  —  -3.394  ft-lbs  and  M2  =  8.345  ft-lbs,  respectively.  All  the  force 
components  increased  in  magnitude  compared  to  the  inviscid  prediction.  The  increase  in  drag 
forces  /  and/;  are  to  be  expected.  The  increase  in  lift  force/  is  due  to  the  higher  peak  pressure 
on  the  bottom  surface  of  1.43  x  10‘3  psi  gauge  (Fig.  13b)  for  the  viscous  case  compared  to  a 
peak  value  of  1.27  x  10‘3  psi  gauge  for  the  inviscid  case  although  the  pressure  distribution  on  the 
top  surface  is  quite  similar  (Fig.  13a).  The  increased  peak  pressure  on  the  bottom  surface  is  due 
to  the  presence  of  the  leading  edge  vortex  that  can  be  seen  in  the  velocity  vectors  on  the  y=0 
plane  shown  in  Fig.  14. 

Effect  Of  The  Presence  Of  Trusses 

The  supporting  trusses  were  added  to  the  reflector  from  the  CAD  model  given  in  terms  of 
IGES  files  from  Code  8230.  These  1GES  descriptions  were  converted  to  fecad  format  using  the 
GridTool6  software.  This  CAD  model  was  modified  so  that  it  is  suitable  for  CFD  calculations  by 
matching  various  members  of  the  trusses  and  eliminating  very  small  gaps  in  the  structure 
(Fig.  15),  A  volume  mesh  was  generated  for  this  configuration  consisting  of  106,310  boundary 
points,  526,636  total  points  and  2,805,589  tetrahedral  elements.  Inviscid  and  viscous  flows  past 
this  configuration  were  computed  for  the  configuration  rotating  at  cn  =  60  rpm  in  the  clockwise 
direction  about  the  y-axis, 

The  convergence  history  of  the  forces  is  shown  in  Fig,  16.  It  is  clear  that  the  differences 
between  inviscid  and  viscous  cases  are  negligible.  The  converged  forces  and  moments  about  the 
origin  shown  in  Fig.  15,  for  the  inviscid  case  are/,  =  -0.303  lbs,  /  =  1.189  lbs  and /  =  0.56  lbs, 
Mx  =  -2.66  ft-lbs,  My  =  -3.50  ft-lbs  and  M,  -  6,30  ft-lbs,  respectively;  for  the  viscous  case,  the 
forces  and  moments  are  /  =  -0.302  lbs,  /  -  1.084  lbs  and  /  =  0.641  lbs,  Mx  =  -2.69  ft-lbs,  My  = 
-3.84  ft-lbs  and  Mz  =  5.74  ft-lbs,  respectively.  Comparing  these  forces  for  the  case  without  the 
trusses,  we  find  that  the  drag  force  in  the  x-direction  (fx )  is  reduced  in  magnitude,  from  -0.99  lbs 
to  ~0.31bs,  due  to  the  influence  of  the  keel  structure  above  the  reflector.  The  drag  force  in  the  z- 
direction  (f. )  is  nearly  unchanged.  The  lift  force  (fy  )  on  the  other  hand  is  reduced  considerably 
from  3.95  lbs  to  1.08  lbs.  To  understand  the  reason  for  this  reduced  lift,  the  pressure 
distributions  on  the  reflector  surfaces  for  the  inviscid  and  viscous  cases  are  plotted  and  are 
shown  in  Figs.  17  and  18  respectively.  Comparing  this  pressure  distribution  to  the  case  without 
the  trusses  (Fig.  9  and  Fig.  13),  we  see  that  the  peak  pressure  on  the  bottom  surface  in  the  case 
without  the  trusses  is  much  higher  and  extends  over  a  larger  area.  Also,  the  region  of  low 
pressure  (blue)  on  the  bottom  surface  in  the  case  with  the  trusses  is  much  larger  leading  to  a 
lower  lift  force.  The  influence  of  the  trusses  on  the  pressure  can  be  clearly  seen  on  the  top 
surface. 

The  moments  about  the  origin  shown  in  Fig.  15  are  computed  for  the  case  of  the  standalone 
reflector.  The  local  origin  is  displaced  from  the  origin  of  the  entire  system  (reflector  and  trusses) 
mainly  in  the  y  direction.  The  displacements  in  the  x  and  y  directions  are  0.86"  and  82,9", 
respectively.  Hence,  the  moments  about  the  origin  of  the  entire  system  can  be  computed.  For  the 
inviscid  case,  they  are  Mx  =  -2.001  ft-lbs,  My  -  -3.226  ft-lbs  and  Mt  =  15.231  ft-lbs,  and  for  the 
viscous  case,  Mx  =  -1.612  ft-lbs,  My  =  -3.441  ft-lbs  and  M,  =  17,286  ft-lbs.  The  moments  about 
the  x-  and  y-axes  are  of  the  same  order  of  magnitude  in  comparison  with  the  trusses,  but  the 
moments  about  the  2-axis  differs  considerably.  This  change  is  due  to  the  partly  due  to  the  change 
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in  the  fx  and  fy  components  of  the  force.  The  other  contribution  to  the  change  in  Mz  arises  due  to 
the  change  in  center  of  pressure.  The  center  of  pressure  for  the  reflector  was  computed  for  the 
viscous  cases  with  and  without  the  trusses.  For  the  case  without  the  trusses,  the  center  of 
pressure  was  located  at  xp=  (30.01",  70.79",  16.82").  In  the  case  with  trusses,  the  center  of 
pressure  was  found  to  be  at  xp  =  (50.76",  45.83",  56.91").  The  change  in  the  x  and  y  moment 
arms  thus  contribute  to  the  majority  of  the  change  in  the  moment  component  Mv  The  forces  and 
moments  for  all  the  cases  computed  were  provided  to  the  Spacecraft  Engineering  Department  for 
computations  of  the  jitter  angle  of  the  rotating  system. 

The  flow  between  the  reflector  and  the  base  of  the  supporting  trusses  is  shown  in  Fig.  19. 
Particle  traces  near  the  trailing  edge  of  the  reflector  are  shown  in  Fig.  19a.  These  particles  are 
released  from  a  rake  of  rectangular  grid  of  points  (shown  in  light  green  color  in  Fig.  19a)  whose 
plane  is  parallel  to  the  diameter  of  the  reflector  and  0.5"  displaced  from  the  edge  of  the  reflector. 
Using  the  instantaneous  velocity  field,  the  positions  of  these  particles  were  obtained  by 
integrating  the  velocity  at  these  points  until  the  length  of  these  traces  exceeded  a  specified 
length,  or  the  particles  ended  on  a  solid  boundary,  or  exited  the  computational  domain.  These 
particles  are  colored  according  to  the  magnitude  of  velocity  (in  inches/sec.)  at  that  location.  A 
trailing  edge  vortex  can  be  seen  from  the  particle  traces  in  Fig.  19a.  Die  flow  below  the  reflector 
and  the  base  of  the  supporting  truss  consists  of  a  region  of  low  velocity  fluid.  Particles  released 
from  just  below  the  reflector,  Fig.  19b,  show  that  the  values  of  velocity  are  in  the  range  of  0  to 
80  in  s'5.  Also,  these  particles  end  up  on  the  recirculation  region  located  at  the  bottom  of  the 
reflector.  The  contours  of  the  magnitude  of  velocity  on  the  y=0  plane  also  shows  that  the  flow  in 
the  region  between  the  reflector  and  the  base  is  nearly  stagnant. 

ACKNOWLEDGEMENTS 

This  work  was  supported  by  Code  8200  of  NRL  under  a  contract  from  SPAWAR  on  the 
WindSat  project.  The  author  would  like  to  thank  Mr.  Michael  Mook  and  Mr.  Dave  Spencer  of 
Code  8200,  for  their  support  and  Mr.  John  Montgomery  and  Mr.  Andy  Feerst  of  Code  8200  for 
their  help  in  providing  the  geometry  of  the  reflector  and  trusses.  The  valuable  discussions  with 
and  Mr.  John  Gardner  of  Code  6440  and  Prof.  Rainald  Lohner  of  George  Mason  University 
throughout  the  course  of  this  project  is  also  greatly  appreciated.  This  work  was  supported  in  part 
by  a  grant  of  HP  C  time  from  the  DoD  HPC  centers,  ARL  MSRC  SG1-02K  and  NRL  SGI-02K. 


Computations  of  Forces  on  WirtdSat  Reflector 


7 


REFERENCES 

1.  Ramamurti,  R.,  Sandberg,  W.C.  and  Lohner,  R.,  “Simulation  of  a  Torpedo  Launch  Using  a  3- 
D  Incompressible  Finite  Element  Flow  Solver  and  Adaptive  Remeshing,”  AIAA  Paper  No. 
95-0086,  January  1995. 

2.  Ramamurti,  R.  and  Lohner,  R.,  “Evaluation  of  an  Incompressible  Flow  Solver  Based  on 
Simple  Elements,”  Advances  in  Finite  Element  Analysis  in  Fluid  Dynamics,  1992,  FED 
Vol.  137,  Editors:  M.N  Dhaubhadel  et  al. ,  ASME  Publication,  New  York,  pp.  33-42. 

3.  Ramamurti,  R,  Lohner,  R.  and  Sandberg,  W.C.,  “Evaluation  of  Scalable  3-D  Incompressible 
Finite  Element  Solver,”  AIAA  Paper  No.  94-0756, 1994. 

4.  Schlichting,  H.,  Boundary  Layer  Theory,  Seventh  Edition,  McGraw-Hill,  New  York. 

5.  Smagorinsky,  J.,  General  Circulation  Experiments  with  the  Primitive  Equations  I.  The  Basic 
Experiment;  Mon.  Weather  Rev.,  91,  pp.  99-164,  1963. 

6.  Samareh,  J.,  “GridTool:  A  surface  Modeling  and  Grid  Generation  Tool,”  Proc.  of  the 
workshop  on  Surface  Modeling,  Grid  Generation,  and  Related  Issues  in  CFD  Solutions, 
NASA  CP-3291,  May  1995. 


8 


Ravi  Ramamurti 


Fig.  1 .  Flow  near  a  rotating  disk,  from  H.  Schlichtrng. 


Fig.  2.  Computational  domain  for  flow  near  a  rotating  disk, 
a>=60  rpm,  r,  =  1 .0,  r2  =  5.0  and  z^,  =  5. 


Fig.  3.  Comparison  of  velocity  components  with  the  exact  solution. 


Computations  of  Forces  on  WindSat  Reflector 


9 


a.  radial  component 


c.  axial  component 


Fig.  4.  Velocity  and  pressure  distribution  ony=0  plane  for  flow  near  a  rotating  disk,  to=60  ipm. 


Fig.  5.  Flow  between  two  concentric  rotating  cylinders,  inner  radius  r,  -  1,  outer  radius  t2  =3. 
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Fig.  6.  Comparison  of  circumferential  velocity  distribution. 


a.  circumferential  velocity 


b.  pressure  c,  magnitude  of  velocity 


Fig.  7.  Flow  between  two  concentric  rotating  cylinders,  <n=60  rpm. 
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Fig.  8.  Computational  domain  of  the  WindSat  reflector. 
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Fig.  9.  Surface  pressure  distribution  on  the  reflector,  co  =  60  ipm. 
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a.  top  surface 


a.  bottom  surface  b.  close-up  view 


Fig.  11.  Velocity  vectors  on  the  bottom  surface  of  the  reflector,  to  =  60  rpm. 
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Fig.  12.  Convergence  history  of  force  components  for  inviscid  and  viscous  flows. 


a.  top  surface 
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Fig.  13.  Pressure  distribution  on  the  reflector.  Re  =  1.2  x  10s. 
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Fig.  14.  Velocity  vectors  on  y  =  0  plane,  Re  =  1.2  x  10e. 


Fig.  15.  WindSat  reflector  with  supporting  trusses. 
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a.  front  view  b.  top  view 

Fig.  17.  Pressure  distribution  over  the  WindSat  rotating  at  60  rpm  (inviscid). 


a.  front  view  b.  top  view 

Fig.  18.  Pressure  distribution  over  the  WindSat  rotating  at  60  rpm  (viscous). 
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b,  particle  traces  below  the  reflector 


c.  magnitude  of  velocity  on  y=0  plane 
Fig*  19.  Flow  between  the  reflector  and  trasses. 


